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The annihilations of WIMPs produce high energy gamma-rays in the final 
state. These high energy gamma-rays may be detected by imaging atmospheric 
Cherenkov telescopes (IACTs). Amongst the plausible targets are the Galactic 
Center, the centre of galaxy clusters, dwarf Sphreroidal galaxies and substruc- 
tures in Galactic haloes. I will review on the recent results from observations 
of ongoing IACTs. 
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1. Introduction 

Cosmological and astrophysical probes suggest that ^23% of the Universe 
is composed of non-baryonic dark matter (DM), commonly assumed to be 
in the form of Weakly Interacting Massive Particles (WIMPs) arising in ex- 
tensions of the Standard Model of Particle Physics (for reviews see, e.g. 1,2 ). 
Amongst the most widely discussed DM candidates are the lightest neu- 
tralino in supersymmetric extensions of the Standard Model 3 and the first 
excitation of the Kaluza-Klein bosons (LKP) in universal extra dimension 
theories. 4-6 

The annihilation of WIMP pairs can produce in the final state a con- 
tinuum of gamma-rays whose flux extends up to the DM particle mass, 
from the hadronization and decay of the cascading annihilation products. 
In supersymmetric models, the gamma-ray spectrum from neutralino anni- 
hilation is not uniquely determined and the branching ratios (BRs) of the 
open annihilation channels arc not determined since the DM particle field 
content is not known a priori. In contrast, in Kaluza-Klein scenarios where 
the lightest Kaluza-Klein particle (LKP) is the first KK mode of the hyper- 
charge gauge boson, the BRs of the annihilation channels can be computed 
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given that the field content of the DM particle is known. The gamma-ray 
flux from annihilations of DM particles of mass muM accumulating in a 
spherical DM halo can be expressed in the form : 



d<5>(An,E 7 ) _ 1 {av) dN y 
dE-y 8-7T m 2 DM dE-y 



J {An) An (1) 

s v ' Astrophysics 

Particle Physics 

as a product of a particle physics component with an astrophysics compo- 
nent. The particle physics part contains (av), the velocity-weighted anni- 
hilation cross section, and dN-y/dE-y, the differential gamma-ray spectrum 
summed over the whole final states with their corresponding branching 
ratios. The astrophysical part corresponds to the line-of-sight-integrated 
squared density of the DM distribution J, averaged over the instrument 
solid angle An usually matching the angular resolution of the instrument : 

J = ( p 2 {r[s])ds J(AO) = -!- / PSFxJdn (2) 

Ji.o.s An j An 

where PSF stands for the point spread function of the instrument. 

The annihilation rate being proportional to the square of the DM density 
integrated along the line of sight, regions with enhanced DM density are 
primary targets for indirect searches. Among them are the Galactic halo, 
external galaxies, galaxy clusters, substructures in galactic haloes, and the 
Galactic Center. We report here on recent results on dark matter searches 
with current lACTs such as H.E.S.S. and MAGIC, towards the Galactic 
Center, dwarf galaxies from the Local Group, Galactic globular clusters 
and DM mini-spikes around intermediate mass black holes (IMBHs). 



2. The Galactic Centre 

H.E.S.S. observations towards the Galactic Center have revealed a bright 
pointlike gamma-ray source, HESS J1745-290, 7 coincident in position with 
the supermassive black hole Sgr A*, with a size lower than 15 pc. Diffuse 
emission along the Galactic plane has also been detected 8 and correlates 
well with the mass density of molecular clouds from the Central Molecular 
Zone, as traced by CS emission. 9 According to recent detailed studies, 10 the 
source position is located at an angular distance of 7.3"±8.7" tat ±8.5g' yst 
from Sgr A*. The pointing accuracy allows to discard the association of 
the very high energy (VHE) emission with the center of the radio emission 
of the supernova remnant Sgr A East but the association with the pulsar 
wind nebula G359.95-0.04 can not be ruled out. From 2004 data set, the 
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energy spectrum of the source is well fitted in the energy range 160 GeV - 
30 TeV to a power-law spectrum dN/dE oc E _r with a spectral index T = 
2.25±0.04 sta t.±0.1 sys t.. No deviation from a power-law is observed leading 
to an upper limit on the energy cut-off of 9 TeV (95% C.L.). The VHE 
emission from HESS J1745-290 does not show any significant periodicity 
or variability from 10 minutes to 1 year. 10 Besides plausible astrophysical 
origins (see e.g. 11 and references therein), an alternative explanation is the 
annihilation of DM in the central cusp of our Galaxy. The spectrum of HESS 
J1745-290 shows no indication for gamma-ray lines. The observed gamma- 
ray flux may also result from secondaries of DM annihilation. The left 
hand side of Fig. lshows the H.E.S.S. spectrum extending up to masses of 
about 10 TeV, which requires large neutralino masses (>10 TeV). They are 
unnatural in phenomenological MSSM scenarios. The Kaluza-Klein models 
provide harder spectra which still significantly deviate from the measured 
one. Non minimal version of the MSSM may yield flatter spectrum with 
mixed 70% bb and 30% r + r _ final states. Even this scenario does not fit 
to the measured spectrum. The hypothesis that the spectrum measured by 
H.E.S.S. originates only from DM particle annihilations is highly disfavored. 




Fig. 1. Left : Spectral energy density E 2 dN/dE of gamma-rays from HESS J1745- 
290 for 2003 (red empty circles) and 2004 (black filled circles) datasets of the H.E.S.S. 
observation of the Galactic Center. The shaded area shows the best power-law fit to the 
2004 data points. The spectra expected from the annihilation of a MSSM-like 14 TeV 
neutralino (dashed green line), a 5 TeV KK DM particle (dotted pink line) and a 10 TeV 
DM particle annihilating into 70% bb and 30% t+t - in final state (solid blue line) are 
presented. Right: Reconstructed VHE gamma-ray energy spectrum of the GC (statistical 
errors only) as measured by the MAGIC collaboration. The full line shows the result of 
a power-law fit to the data points. The dashed line shows the 2004 result of the HESS 
collaboration. 13 The dot-dashed line shows the energy spectrum of the Crab nebula as 
measured by MAGIC. 
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The MAGIC observations were carried out towards the Galactic Cen- 
tre since 2004 and revealed a strong emission. The observed excess in the 
direction of the GC has a significance of 7.3 standard deviations 12 and is 
compatible with a pointlike source. Large zenith observation angle (>60°) 
implies an energy threshold of ^400 GeV. The source position and the flux 
level are consistent with the measurement of HESS 13 within errors. The 
right hand side of Fig. I shows the reconstructed VHE gamma-ray energy 
spectrum of the GC after the unfolding with the instrumental energy reso- 
lution. The differential flux can be well described by a power law of index 
T = 2.2 ± 0.2 sta t. ± 0.2 syst .. The systematic error is estimated to be 35% 
in the flux level determination. The flux level is steady within errors in the 
time-scales explored within these observations, as well as in the two year 
time-span between the MAGIC and HESS observations. An interpretation 
in term of DM require a minimum value of mass higher than 10 TeV. Most 
probably, if DM signal exists is overcome by other astrophysical emitters. 

3. Dwarf Galaxies from the Local Group 

Dwarf spheroidal galaxies in the Local Group are considered as privileged 
targets for DM searches since they are among the most extreme DM- 
dominatcd environments. Measurements of roughly constant radial velocity 
dispersion of stars imply large mass-to-luminosity ratios. Nearby dwarfs are 
ideal astrophysical probes of the nature of DM as they usually consist of 
a stellar population with no hot or warm gas, no cosmic ray population 
and little dust. Indeed, these systems are expected to have a low intrinsic 
gamma-ray emission. This is in contrast with the Galactic Center where 
disentangling the dominant astrophysical signal from possible more exotic 
one is very challenging. Prior to the year 2000, the number of known satel- 
lites was eleven. With the Sloan Digital Sky Survey (SDSS), a population 
of ultra low-luminosity satellites has been unveiled, which roughly doubled 
the number of known satellites. Among them are Coma Berenices, Ursa Ma- 
jor II and Willman 1. IACTs have started observation campaigns on dwarf 
galaxies for a few years. Table 1 presents a possible list of the preferred 
targets for DM searches. 

The star velocity dispersions in Draco reveal that this object is dom- 
inated by DM on all spatial scales and provide robust bounds on its DM 
profile, which thus decreases uncertainties on the astrophysical contribu- 
tion to the gamma-ray flux. The MAGIC collaboration searched for a 
steady gamma-ray emission from the direction of Draco. 17 The analysis 
energy threshold after cuts is 140 GeV. No significant excess is found. 
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Table 1. A tentative list of preferred dwarf galaxies. 



Name 


Distance 3, 


Luminosity 3 


M/L b 


Best positioned IACTs 




(kpc) 


(10 3 Lq) 


(M /L ) 




Carina 


101 


430 


40 


HESS, CANGAROO 


Coma Berenices 


44 


2.6 


450 


MAGIC, VERITAS 


Draco 


80 


260 


320 


MAGIC, VERITAS 


Fornax 


138 


15500 


10 


HESS, CANGAROO 


Sculptor 


79 


2200 


7 


HESS, CANGAROO 


Sagittarius 


24 


58000* 


25 


HESS, CANGAROO 


Sextans 


86 


500 


90 


HESS, CANGAROO 


Ursa Minor 


66 


290 


580 


MAGIC, VERITAS 


Ursa Major II 


32 


2.8 


1100 


MAGIC, VERITAS 


Willman 1 


38 


0.9 


700 


MAGIC, VERITAS 



Note: a See Ref. 14. b See Ref. 15. *Scc Rcf. 16. 



For a power law with spectral index of 1.5, typical for a DM annihi- 
lation spectrum, and assuming a pointlikc source, the 2a upper limit is 
$(E > 140GcV) = 1.1 x 10 _11 cm _2 s _1 . The measured flux upper limit is 
several orders of magnitude larger than predicted for the smooth DM dis- 
tribution in mSUGRA models. The limit on the flux enhancement caused 
by high clumpy substructures or a black hole is around O(10 3 — 10 9 ). The 
WHIPPLE collaboration has also observed Draco. 18 Fig. 3 shows the 95% 
C.L. exclusion curve on av. The upper limit is at the level of 10~ 22 cm 3 s _1 
for 1 TeV neutralino annihilating with BRs of 90% in bb and 10% in t + t~. 

The Sagittarius (Sgr) dwarf galaxy, one of the nearest Galaxy satellites 
of the Local Group, has been observed by H.E.S.S. since 2006. The anni- 
hilation signal from Sgr is expected to come from a region of 1.5 pc, which 
is much smaller than the H.E.S.S. point spread function. Thus, a pointlikc 
signal has been searched for. No significant gamma-ray excess is detected 
at the nominal target position. A 95% C.L. upper limit on the gamma-ray 
flux is derived: $ 95%C L (E 7 > 250GcV) = 3.6 x RT^cm"^- 1 , assuming 
a power-law spectrum of spectral index 2. 2. 19 A modelling of the Sgr DM 
halo has been carried out. Two models of the mass distribution for the 
DM halo have been studied: a cusped NFW profile and a cored isothermal 
profile, to cmcompass a large class of profiles. The cored profile has a small 
core radius due to a cusp in the luminous profile. The value of the line- 
of-sight-integrated squared density is then found to be larger for the cored 
profile than for the NFW profile (see Ref. 19 for more details). The left 
hand side of Fig. 2 presents the constraints on the velocity-weighted anni- 
hilation cross section av for a cusped NFW and cored profiles in the solid 
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angle integration region Af2 = 2 x 10~ 5 sr, for neutralino DM. Predictions 
for SUSY models are displayed. For a cusped NFW profile, H.E.S.S. does 
not set severe constraints on av. For a cored profile, due to a higher central 
density, stronger constraints are derived and some pMSSM models can be 
excluded in the upper part of the scanned region. Although the nature of 
Canis Major is stilll debated, HESS observed this object and results on 
Canis Major 20 are presented in Ref. 21. 
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Fig. 2. Left: Upper limits at 95%C.L. on av versus the neutralino mass for a cusped 
NFW and cored DM halo profile for Sgr. The predictions in pMSSM are also plotted 
with in addition those satisfying in addition the WMAP constraints on the cold DM 
density. Right: Willman 1 a-plot as seen by MAGIC in 15.5 hours above a fiducial energy 
threshold of 100 GeV. The red crosses represent the ON— data sample, the blue shaded 
region is the OFF-data sample normalized to the ON-data sample between 30° — 80° . 
The vertical red dotted line represents the fiducial region a < 12° where the signal is 
expected. 

The recently discovered Willman 1 dwarf galaxy has been observed by 
MAGIC. 22 Willman 1 has a a total mass of -^5xl0 5 M Q , which is about 
an order of magnitude smaller than those of the least massive satellite 
galaxies previously known. This object has one of the highest mass-to- 
luminosity ratio. No significant gamma-ray excess beyond 100 GeV above 
the background was observed in 15.5 hours of observation of the sky region 
around Willman 1. This is shown in the right hand side of Fig. 2, where the 
"a-plot" is reported a . The signal is searched with a cut slightly larger than 
for a pointlike source to take into account a possible source extension which 
makes a fiducial region to be a <12°. Flux upper limits of the order of 10~ 12 
cm _2 s _1 for benchmark mSUGRA models are obtained. 22 Boost factors in 
flux of the order of 10 3 are required in the most optimistic scenarios, i.e. in 



a The o-parameter is the angular distance between the shower image main axis and the 
line connecting the image barycenter and the camera center. 
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the funnel region or in region where the internal bremsstrahlung may play 
an important role. However, uncertainties on the DM distribution or the 
role of substructures may significantly reduce this boost. 



4. Galactic Globular Clusters 

Globular clusters are not believed to be DM-dominated objects. Their mass- 
to-luminosity ratios are well described by purely King profile which suggets 
no significant amount of DM. However, the formation of globular clusters 
fits in the hierarchical strucure formation scenario in which globular clus- 
ters may have formed in DM overdensities. During their evolution, they 
may have hold some DM in their central region. M15 is a nearby galac- 
tic globular cluster with a core radius of 0.2 pc and a central density of 
^10 7 M Q pc~ 3 . Even if the mass-to-light ratio of M15 does not suggest a 
significant component of DM, the compact and dense core of M15 does not 
prevent from looking for a hypothetical DM signal. The presence of stars 
and DM in globular clusters may lead to a DM enhancement in their inner 
part. This is usually treated with the adiabatic contraction model. 23 

The Whipple collaboration obsered M15 in 2003. The modelling of DM 
profile assuming adiabatic contraction of DM from an initial NFW dis- 
tribution leads to an enhancement in the astrophysical factor of ~10 2 to 



Fig. 3. Exclusion limit on av as a function of the neutralino mass m x on M15 (dotted 
line) assuming an annihilation spectrum with 90% in bb and 10% in t+t~ , and the NFW 
profile after adiabatic contraction. The solid line is the exclusion limit for Draco and the 
dashed line for Ursa Minor. 





M x [GeV] 
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10 3 . Fig. 3 shows the exclusion curve on av (dotted line) versus the neu- 
tralino mass assuming a NFW profile after adiabatic contraction. 18 The 
limit reaches ~10 -24 cm 3 s _1 for TeV neutralinos. 

5. Galaxy clusters 

Astrophysical systems such as the VIRGO galaxy cluster 24 have been con- 
sidered as target for DM annihilation searches. The elliptical galaxy M87 
at the center of the VIRGO cluster has been observed by H.E.S.S. It is an 
active galaxy located at 16.3 Mpc which hosts a black hole of 3.2xlO 9 M . 
From 2003 to 2006, H.E.S.S. detected a 13cr signal in 89 hour observation 
time at a location compatible with the nominal position of the nucleus of 
M87. 25 With the angular resolution of H.E.S.S., the VHE source is point- 
like, with an upper limit on the size of 3' at 99% confidence level. Given 
the distance of M87, this corresponds to a size of 13.7 kpc. 




Energy (TeV) 



Fig. 4. Differential energy spectrum of M87 from ~400 GeV up to 10 TeV from 2004 
and 2005 data. Corresponding fits to a power-law function dN/dE E are plotted as solid 
lines. The data exhibit indices of T = 2.62±0.35 for 2004 and T = 2.22±0.15 for 2005. 

Fig. 4 shows the energy spectrum for 2004 and 2005 data. Both data sets 
are well fitted to a power-law spectrum with spectral indices T — 2.62±0.35 
(2004) and T = 2.22±0.25 (2005). The 2005 average flux is found to be 
higher by a factor ~5 to the average flux of 2004. The integrated flux above 
730 GeV shows a yearly variability at the level of 3.2a. 25 Variability at the 
day time scale has been observed in the 2005 data at the level of 4cr. This 
fast variability puts stringent constraints on the size of the VHE emission 
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region. The position centered on M87 nucleus excludes the center of the 
VIRGO cluster and outer radio regions of M87 as the gamma-ray emission 
region. The observed variability at the scale of ^2 days requires a compact 
emission region below 50R S , where K s ^1CP 15 cm is the Schwarzschild 
radius of the M87 supermassive black hole. 25 The short and long term 
temporal variability observed with H.E.S.S. excludes the bulk of the TeV 
gamma-ray signal to be of DM origin. 

6. Galactic Substructures: the Case for IMBHs 

Mini-spikes around Intermediate Mass Black Holes have been recently pro- 
posed as promising targets for indirect dark matter detection. 26 The growth 
of massive black holes inevitably affects the surrounding DM distribution. 
The profile of the final DM overdcnsity, called mini-spike, depends on the 
initial distribution of DM, but also on astrophysical processes such as gravi- 
tational scattering of stars and mergers. Ignoring astrophysical effects, and 
assuming adiabatic growth of the black hole, if one starts from a NFW 
profile, a spike with a power-law index 7/3 is obtained, as relevant for the 
astrophysical formation scenario studied here characterized by black hole 
masses of ~10 5 M . 26 Mini-spikes might be detected as bright pointlikc 
sources by current IACTs. 

H.E.S.S. data collected between 2004 and 2007 during the Galactic plane 
survey have allowed to accurately map the Galactic plane between ±3° in 
galactic latitude and from -30° to 60° in galactic longitude with respect 
to the Galactic Center position. The study of the H.E.S.S. sensitivity in a 
large field of view to dark matter annihilations has been performed. 2 Fig. 
5 shows the experimentally observed sensitivity map in the Galactic plane 
from Galactic longitudes l=-30° to l=+60° and Galactic latitudes b=-3° 
to b=+3°, for a DM particle of 500 GeV mass annihilating into the 100% 
BR bb channel. The H.E.S.S. sensitivity depends strongly on the exposure 
time and acceptance maps which are related to the choice of the pointing 
positions. The flux sensitivity varies along the latitude and longitude due 
to inhomogeneous coverage of the Galactic plane. In the band between -2° 
and 2° in Galactic latitude, a DM annihilation flux sensitivity at the level 
of 10~ 12 cm _2 s _1 is achieved for a 500 GeV DM particle annihilating in 
the bb channel. Deeper observations of the Galactic Center and at Galactic 
longitude of ^20° allow the flux sensitivity to be of ^5xl0~ 13 cm _2 s _1 . 
For b>2°, the sensitivity is deteriorated due to a weaker effective exposure. 
For b=0° and l=-0.5°, the flux sensitivity is ~10~ 13 cm^s -1 in bb. 

H.E.S.S. reached the required sensitivity to be able to test DM an- 
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Fig. 5. H.E.S.S. sensitivity map in Galactic coordinates, i.e. 90% C.L. limit on the 
integrated gamma-ray flux above 100 GeV, for dark matter annihilation assuming a DM 
particle of mass m x = 500 GeV and annihilation into the bb channel. The flux sensitivity 
is correlated to the exposure and acceptance maps. In the Galactic latitude band between 
-2° and 2°, the gamma-ray flux sensitivity reaches 10 — 12 cm~ 2 s _1 . 



nihilations from mini-spikes in the context of one relatively favorable sce- 
nario for IMBH formation and adiabatic growth of the DM halo around the 
black hole (e.g. scenario B of Ref . 26). H.E.S.S. observations (2004-2006) 
of the Galactic plane allowed to discover more than 20 very high energy 
gamma-ray sources . 28 Some of them have been identified owing to their 
counterparts at other wavelengths, but almost half of the sources have no 
obvious counterpart and are still unidentified. 29 An accurate reconstruc- 
tion of their energy spectra shows that all the spectra are consistent with 
a pure power-law, spanning up to two orders of magnitude in energy above 
the energy threshold. None of them exhibits an energy cutoff, characteris- 
tic of DM annihilation spectra, in the energy range from ^100 GeV up to 
10 TeV. Furthermore, the detailed study of their morphology 29 shows that 
all the sources have an intrinsic spatial extension greater than ~5 arcmin- 
utes, while mini-spikes are expected to be pointlike sources for H.E.S.S. No 
IMBH candidate has been detected so far by H.E.S.S. within the survey 
range. Based on the absence of plausible IMBH candidates in the H.E.S.S. 
data, constraints are derived on the scenario B of Ref. 26 for neutralino 
or LKP annihilations, shown as upper limits on av. 27 Fig. 6 shows the ex- 
clusion limit at the 90% C.L. on av as a function of the neutralino mass. 
The neutralino is assumed to annihilate into bb and t + t~ with 100% BR, 
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respectively. Predictions for SUSY models are also displayed. The limits on 



*rio _ 

E 

■ — <it\-ii - 



— i — i — i — i — i 1 1 1 1 



10 



— i — i — i — i — i 1 1 1 

H.E.S.S.90%C.L. - 
pMSSM 

pMSSM + WMAP 




m DM [TeV] 



Fig. 6. Constraints on the IMBH gamma-ray production scenario for different neu- 
tralino parameters, shown as upper limits on av as a function of the mass of the neu- 
tralino mrjAfi brrt with a number of implicit assumptions about the IMBH initial mass 
function and halo profile (see Ref. 27 for details). For the scenario studied here, the 
probability of having no observable haloes in our Galaxy is 10% from Poisson statistics 
making these limits essentially 90% C. L. exclusion limits for this one particular (albeit 
optimistic) scenario (grey shaded area). The DM particle is assumed to be a neutralino 
annihilating into bb pairs or r + r — pairs to encompass the softest and hardest annihila- 
tion spectra. The limit is derived from the H.E.S.S. flux sensitivity in the Galactic plane 
survey within the mini-spike scenario. SUSY models (black points) arc plotted together 
with those satisfying the WMAP constraints on the DM particle relic density (magenta 
points). 

(TV are at the level of 10 -28 cm 3 s _1 for the bb channel for neutralino masses 
in the TeV energy range. Limits are obtained one mini-spike scenario and 
constrain on the entire gamma-ray production scenario. 



7. Conclusion 

Dark matter searches will continue and searches with the phases 2 of 
H.E.S.S. and MAGIC will start soon. The phase 2 of H.E.S.S. will consist 
of a new large 28 m diameter telescope located at the center of the existing 
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array, and will allow to lower the analysis energy threshold down to less 
than 50 GeV. The installation of the second telescope on the MAGIC site 
will allow to perform more sensitive searches by the use of the stereoscopic 
mode. The upcoming generation of IACTs, the Cherenkov Telescope Ar- 
ray, is in the design phase. This array composed of several tens of telescopes 
will permit to significantly improve the performances on both targeted DM 
searches and widc-field-of-vicw DM searches. 
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